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Fundamental physics with neutron stars 
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Abstract Neutron stars are rich laboratories of multiple branches of modern physics. These include gravi- 
tational physics, nuclear and particle physics, (quantum) electrodynamics, and plasma astrophysics. In this 
chapter, we present the pioneering theoretical studies and the pivotal historical observations on which our 
understanding of neutron stars is based on. Then, we discuss the usage of neutron stars as probes of funda- 
mental theories of physics. 
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Introduction 


Neutron stars are extremely dense, rapidly rotating, superfluid magnets. They shine to us across a broad 
electromagnetic spectrum from low-frequency radio bands to high-energy X-rays and gamma-rays. Some 
of them produce continuous radiation, some show pulsations, and some explode in flares. This activity is 
powered by their huge gravitational, magnetic, rotational, and chemical energy. These enormous reservoirs 
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of energy, in turn, drive exotic phenomena encompassing many fields of modern physics, from the quantum 
realm to space-time-deforming general relativity. Most of the observed phenomena from neutron stars are 
so extreme that they can not be studied in contemporary terrestrial laboratories. This makes neutron stars the 
cosmic laboratories of extreme physics. 

A canonical neutron star has a radius of R ~ 12km and a mass of M ~ 1.5Mo; the observed masses 
range from about 1.1 to 2.1 Mo, and theoretically supported radii range from about 10 to 15km. This makes 
them extremely dense astrophysical objects capable of significantly distorting the local space-time—hence, 
general relativity is needed when describing the physical laws around them. Neutron stars are also observed 
to rotate with spin periods ranging from P ~ 0.001 to 10s. When such a tightly compressed object is rapidly 
rotating, it stores massive amounts of angular momentum that will result in additional “twisting” of the 
space-time “fabric”, in addition to the more usual curving effect caused by the concentrated mass. The first 
section of this chapter describes the gravitational physics of neutron stars. 

The large density means that the matter inside neutron stars is under tremendous pressure. The mean 
density inside the star, Ø ~ 10!5 gcm~?, is comparable to the density inside the atomic nuclei. The matter is 
compressed into such a small volume that quantum mechanical effects define its behavior. How the matter 
responds to this squeezing, as the particles try to counter the compressing force, is encapsulated in the—still 
unknown—equation of state (EOS) of the cold, dense matter. The nuclear and particle physics of neutron 
stars are described in the second section of this chapter. 

Neutron stars have strong magnetic fields, ranging from B ~ 10° to 10!°G. As the star spins, the topology 
of this dynamic magnetosphere changes. The physics of the magnetosphere is governed by both electromag- 
netism, describing the dynamical evolution of the electric and magnetic fields, and quantum electrodynamics 
(QED), describing the coupling of the electromagnetic fields and particles. The electrodynamics, QED phe- 
nomena, and plasma physics of neutron stars are discussed in the third and fourth sections of this chapter. 


Formation of neutron stars 


The idea of stars composed of neutrons originates from Lev Landau [80, 163]; he was the first to speculate 
on the possibility of squeezing the matter inside regular stars so tightly that quantum effects would start to 
dominate: ”...the density of the matter becomes so great that atomic nuclei come in close contact, forming 
one gigantic nucleus”. Interestingly, this idea was coined even before the existence of neutrons was known.” 
Similar ideas were later proposed also by Walter Baade and Fritz Zwicky [1 2]—they were even bold enough 
to suggest that neutron stars might form in supernovae. This turned out to be exactly right! 

Neutron stars are now known to be born from a collapse of a main-sequence star, slightly more massive 
than our own Sun, with a mass M, ~ 10Mọ, where Mọ © 1.99 x 10°? g [141]. Such a star has a radius 
of approximately 5x that of our own Sun, R, ~ 5Rọ where Ro œ~ 6.96 x 10!°cm.? Then, the resulting 
(very approximate) estimate for the mean density of the matter inside the progenitor star is p, + 0.1 gcm™?, 
roughly b of the density of liquid water. 

Such stars fuse lighter elements into heavier metals, all the way up to iron. Eventually, the heavy iron 
elements accumulate and sink to the bottom of the gravitational potential, forming an iron core for the star. 
This core has a theoretically known maximum mass of about 1.4Mọo. Once this limit is exceeded, the core 
collapses, and the outer layers explode as a bright supernova. 


? Existence of electrically neutral particles with a mass my ~ Mp 5X 10~* g was proven by Chadwick in 1932 [32, 33]. 


3 In reality, the massive star evolves off the main sequence before the supernova explosion and can expand in its radius. However, 
the discussed values are good “first guesses” to perform the order-of-magnitude estimates presented in this section. 
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The remnant core consists of roughly M /mp ~ Ma /mp ~ 10°” atoms, where Mp © 1.67 x 10-*4 ¢ is the 
mass of a proton. Suppose those atoms shed their electrons and form a giant blob of neutrons. In that case, 
the resulting remnant will have a characteristic size 0 ~ (1057) 37, ~ 10km, where ra ~ 107 cm is the 
approximate size of an atomic nucleon. A compact object with a mass ~ Mo and radius ~ 10km is hence 
left behind— a neutron star! Even though these values are approximate and based on very naive conservation 
laws, they are not far from the actual observed values for neutron stars; here we adopt the canonical values 
of a neutron star mass Mys = 1.4Mo and radius Rys = 12km. 

A canonical neutron star has an extraordinarily high mean density, ð ~ Mys [Ris ~ 10 gem~? and 
surface gravity, g ~ 10'4cms~?. The central density of such an object, Pe ~ Pp, far exceeds the typical 
density of a nucleon, Py ~ mp/ r? ~ 10!4gcm~3 < pe. This implies that the matter is compressed so densely 
that there is no empty space between the atomic core and the electrons; instead, individual atomic nuclei are 
packed side by side, forming a macroscopic “quantum ball” of neutrons. 


3 


First observation of a neutron star 


While the theoretical existence of neutron stars was speculated already in the 1920s, conclusive observational 
evidence came only at the end of the 1960s. Neutron stars were first detected using a radio antenna in 
the pastures of Cambridge [63]. The primitive antennas were originally designed for monitoring the radio 
emission from quasars; they were, however, sensitive enough to detect all kinds of human-made noise signals 
(arc welders, sparking thermostats, etc.) Among many such noise signals, Jocelyn Bell, a student of Tony 
Hewis, noticed a signal that did not seem to have a direct connection to anything happening on Earth; instead, 
the signal — that was sometimes on, sometimes off — seemed to originate from a specific location in the sky. 
The mystery deepened when they were able to resolve the signal, using a high-speed recorder, into a regular 
set of pulses with a stable period of P ~ 1.337s. The source was initially (and jokingly?) named LGM-1, 
standing for Little Green Men, since they did not know of any astrophysical engine capable of powering this 
periodic extraterrestrial signal. 

By 1968, the radio observers had detected multiple pulsing signals with periods varying from P ~ 0.25 to 
3.3s. The sub-structure of the pulses themselves could be even shorter, At < 0.01s. The sources were also 
found to be extremely stable. There was, however, a small systematic decrease detected in many of them, 
on a timescale of P/P ~ 10° to 10° yr, where P = dP/dt is the so-called period-derivative, quantifying the 
change of the period.* 


Theoretical arguments for the existence of neutron stars 


So what can we conclude from these observational facts? The following reasoning was laid out by an Amer- 
ican astrophysicist Thomas Gold [52]. 

Firstly, the short duration of the pulses implies a tiny object; light crosses a distance of ~ 3000km in 
a duration of the sub-pulse, At ~ 0.01s. This means that the emitting region must be smaller than that. In 


4 Such highly accurate measurements were possible because the signals were very stable and could be accumulated over long 
time intervals of many years. 


4 Joonas Nattila and Jari J.E. Kajava 


practice, only black holes, neutron stars, or white dwarfs are small enough astrophysical objects so that their 
radius R < 3000km. This was the first clue: the culprit has to be some type of compact object. 

Secondly, how to produce a stable pulsation from such sources? The usual suspects that could introduce 
such highly periodic signals include rotation, vibration, or orbital motions. Orbital motions are, however, 
ruled out because the observed periods of the signal were not increasing (corresponding to inspiral) but were 
found to be decreasing (i.e., negative period derivatives). In addition, a binary system with an orbital period 
of Porb ~ 1s would merge in a few years as it loses angular momentum via gravitational radiation. Therefore, 
orbital motions could not be the source of the pulsations. 

The shortest time scale out of the proposed emission mechanisms results from vibrations; here, the gravity 
of the object is a natural candidate for the restoring force. The timescale for this is £ ~ 2% //Gp, which, for 
the densities of white dwarfs, p ~ 108 gcm~?, is ~ 1s; borderline too slow to explain the fastest observed 
pulse periods. Conversely, for neutron stars, corresponding to characteristic densities p ~ 10'4 gem~?, the 
vibration timescale would be milliseconds; borderline too fast to explain the slowest observed pulses. Black 
holes do not have a surface capable of emitting vibrations, so they do not need to be considered for this 
mechanism. These conflicting timescales were used to argue against vibration as a source of the pulsations. 
Furthermore, it is also quite unlikely that any surface vibration could go on for years and years without 
damping. 

This left the contemporary researchers only with rotation. A maximum rotation rate that a star can support 
before being ruptured will be close to its Keplerian frequency, fx = (1/27) GM /R3; it is fx ~ 1700 Hz for 
a neutron star and fx ~ 0.1 Hz for a white dwarf (with Mwp ~ Mo and Rwp ~7 x 108 cm). This, in the end, 
helped rule out white dwarfs from the possible list of sources because they would have a hard time sticking 
together under such a rapid rotation. Any surface emission from rotating black holes was also naturally ruled 
out because black holes do not have a surface. 

This reasoning led contemporary scientists to believe that the Little Green Men observed in the Cambridge 
pasture originated from rotating neutron stars. 


Rotating neutron stars 


So, what rotation rates can we expect neutron stars to have? The magnitude of the spin can be estimated 
from the conservation of angular momentum during the formation of the proto-neutron star. Our own Sun 
rotates around itself with a period of Ps ~ 25.5 days or Vo œ 4.5 x 1077 Hz. If such a star is compressed 
into a neutron star, and the angular momentum is conserved during the process,> the new spin would be 
v ~ (Ro/Rys)? Vo ~ 2000Hz. This value is a bit high, but it does get us to the right ballpark. Observed spin 
frequencies of neutron stars vary from ~ 0.01 Hz up to ~ 700 Hz [28]. Some X-ray pulsars have been argued 
to rotate as slowly as ~ 1075 Hz [142]. 

The spin period of neutron stars is not constant; even the first pioneering observations of pulsars could 
detect a slow downward drift in the spin period. This spin-down of pulsars is slow and continuous with 
P/P ~ 107!?ss7!, and it results from magnetic dipole losses as the magnetosphere rotates and generates 
electromagnetic radiation which, in turn, carries a fraction of the star’s angular momentum away. This kind 


> A realistic formation of a protoneutron star in a supernova collapse might not conserve angular momentum; for example, 
suppose there is an off-axis stream of material falling toward the center. In that case, it can result in a seeming over-accumulation 
of net angular momentum in the core parts. Hence, estimates of the initial angular momentum of a proto-neutron star are very 
uncertain. 
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of spin-down is, therefore, naturally expected for isolated pulsars. In addition, many pulsars show abrupt 
spin-down events, called glitches, where the spin suddenly changes in seconds. 

There are also observations of neutron stars in binary systems - in fact, a big fraction of stars in any galaxy 
is expected to reside in systems with multiple stars; exotic multiples up to seven stars are known [158]. These 
binary systems are evolving, and the stars commonly exchange mass. When the mass is transferred onto the 
neutron star, it can be spun up: as the in-falling material streams in, hits the star off center, and exerts a 
torque, it can transfer some of the binary system’s angular momentum into the central object. Alternatively, 
the torque can be transferred to the star via the magnetic field lines that are attached to the inner disc. 
This angular momentum transfer generates a sub-population of neutron stars called “recycled millisecond 
pulsars” - neutron stars spun up by accretion to spin periods of up to P ~ 0.001s. 

This also highlights the importance of understanding the surroundings of neutron stars when trying to 
explain their physics. In many cases, the evolution and physics of the neutron stars in binary systems are 
intimately influenced by the companion star. 


Magnetic fields of neutron stars 


The remaining fundamental quantity of a neutron star is its magnetic field. The neutron star’s magnetic 
field strength can be estimated from the conservation of magnetic flux during its birth [160]. A typical 
progenitor of a neutron star has a magnetic field of Bg ~ 1G; for comparison, a typical fridge magnet has 
a magnetic field of ~ 10G. The stellar magnetic fields are powered by dynamo processes in which the 
differential rotation and convective circulations inside the star amplify the magnetic field. The field can be 
highly tangled and irregular due to the turbulent motions of the stellar interiors. 

When such a star collapses into a neutron star, the dynamo process is thought to halt and the internal 
magnetic field to be squeezed into a smaller volume. The resulting strength of the new field can be estimated 
from the conservation of the magnetic flux through the surface (g = fB- dS, where dS is an infinitesimal 
area element); alternatively, we can think that the number of magnetic field lines penetrating the star’s surface 
is conserved. Thus, by conservation of magnetic flux, we end up with Bys ~ Bs (Ro/Rns)* ~ 10!°G. The 
resulting field is enormous; the largest continuously-operating human-made magnets on Earth can sustain 
“only” ~ 10°G. Observed magnetic fields of regular neutron stars range from 10° to 10!*G. 

Some neutron stars are observed to house an even stronger magnetic field of ~ 10!4 to 10!6 G—these 
are known as magnetars [155]. The atypical magnetic field is thought to originate from a dynamo process 
that occurs during the progenitor star’s collapse, amplifying the star’s regular magnetic field by an unknown 
factor of ~ 100 to 1000 to super-high strengths. Additionally, these magnetars are thought to be so young, 
t < 10° yr, that the enormous field has not yet had time to (resistively) decay away. These magnetars were 
discovered because of the most violent astrophysical events ever observed by humankind. 


Gamma-ray blasts from the past 


At 10:51 am EST on the 5th of March, 1979, two unmanned Soviet space probes, Venera 11 and Venera 12, 
suddenly became overloaded with an unprecedented gamma-ray shower.° The gamma rays first hit Venera 


6 See https: //solomon.as.utexas.edu/magnetar.html for more details. 
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11 and then, 5 seconds later, Venera 12. The onboard detectors observed a sudden jump from a regular 
background gamma-ray count rates of ~ 100 up to ~ 40000 counts per second, and then a level entirely off 
the scale, all during a At ~ ms. After this, 11 seconds later, the gamma rays reached an American space 
probe Helios 2 orbiting the Sun, also knocking it out. None of the detectors were designed to withstand such 
intense emission, which saturated the internal photon counters. 

On that day, an unprecedented gamma-ray wavefront continued to sweep through our solar system. After 
the Helios 2 probe, it reached the planet Venus, where the Pioneer Venus spacecraft detected a burst of 
gamma rays and quickly saturated. After seven more seconds, the wavefront reached Earth and was detected 
by Vela satellite and Soviet Prognoz 7, again knocking their internal detectors out. Even the Einstein X-ray 
observatory, orbiting Earth, detected a strong signal, even though it was not pointing in the direction of the 
blast — the gamma rays were diffusing through its back metal shield and still overloading the detector. A 
similar thing happened to the international Sun-Earth Explorer (ISEE) that was pointing in another direction, 
but the gamma rays penetrating through its solid metal body were still energetic enough to ramp up the 
detector to the maximum. 

Fourteen and half hours later, another fainter burst arrived from the same spot, lasting 1.5s. Then, a 
month later, on April 5th and again on April 24th, two more gamma-ray bursts were detected. Humankind 
was registering its first giant flares from what became known as soft gamma-ray repeaters. 

Tracking of the 5th of March wavefront that swept through our Solar system enabled contemporary re- 
searchers to pinpoint the origin of the blast to a nearby supernova remnant, SNR N49. However, this su- 
pernova remnant was not even in our Galaxy but ~ 50kpc away in one of our satellite galaxies, Large 
Magellanic Cloud. The intense blast had traveled almost 160000 years to us and was still powerful enough 
to knock out our space-borne detectors. The nature of the soft gamma-ray repeaters and these giant flares 
remained a mystery for years as no plausible energy source was established to explain the production of the 
immense gamma rays observed. 

Many years later, Robert Duncan and Chris Thompson, both at Princeton University, started investigating 
the implications of hot and messy proto-neutron stars right after they were born in supernova explosions 
[45]. Recent simulations seemed to indicate that the dense fluid in the core had strong convective motions 
that helped to carry the heat out from the hot remnant. The ultra-dense matter at the interior would also 
carry electric currents; hence, a dynamo process similar to Earth’s core could, in theory, be possible. If the 
neutron star would collapse fast enough, rotate rapidly, and have strong convective motions, it could house a 
magnetic dynamo inside of it, which could amplify any existing magnetic field by twisting and folding it to 
anywhere from 10 to 1000 times of its original strength. Therefore, if typical pulsars with B ~ 10!*G were 
the outcome of a failed dynamo, what would happen if the dynamo succeeded? They hypothesized that a 
successful dynamo process would result in a neutron star with an enormous magnetic field up to B ~ 10!°G 
field. These monsters they, later on, dubbed “magnetars”. 

The first observational implication of such a field would be that if the star’s internal stresses twisted and 
shook the magnetosphere, it could exhibit spectacular flares similar to our own Sun. Therefore, in the case 
of a magnetar, the flare could turn a fraction of that immense magnetic field into radiation, and, in theory, 
power intense blasts of gamma-rays with luminosities L > 10® erg s~!, capable of explaining the observed 
giant flares saturating our X-ray and gamma-ray detectors. Eventually, this became the accepted explanation 
for the giant flares, rendering magnetars the kings of the magnetic fields in the known Universe. 
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Many observational faces of neutron stars 


As we have seen, neutron stars have multiple observational ’ faces”. We have detected both isolated neutron 
stars and neutron stars in binary systems. The types of companion stars in the binary systems vary from low- 
mass stars to other neutron stars (or even black holes). Some neutron stars are powered by rotation (pulsars), 
and some only by their internal heat (e.g., central compact objects found in supernova remnants). 

There is a distinct dichotomy in the observed spin distribution: regular pulsars, identified typically as 
relatively young neutron stars, are slowly spinning (P ~ 1s), whereas older neutron stars in binary systems 
have been spun up by accretion, giving rise to what we know as recycled millisecond pulsars. These recycled 
millisecond pulsars can be further divided into rotation-powered (i.e., radio millisecond pulsars, thought 
to operate similarly to regular pulsars but to spin faster) and accretion-powered (i.e., pulsating on X-ray 
frequencies due to the in-falling material that is channeled to the magnetic poles, creating hot spots onto the 
surface) sources. A few stars, known as transitional millisecond pulsars, are observed to transition between 
these two states. 

Some of the rotation-powered millisecond pulsars are also identified as “spiders” (named so because 
these stars are acting nasty on their companion stars). Spiders generally have irregular radio eclipses due to 
intra-binary material driven off the donor by a pulsar wind. They can be further divided into black widows 
(characterized by very low-mass companions of M, < 0.06M,) and redbacks (characterized by hydrogen- 
rich companions of Me Z 0.1Mo). 


~N 


Various classes of neutron stars and their observational phenomena are listed in Table 1. 


Table 1 Different observational classes of neutron stars 


Class Description # Observed phenomena 
Isolated neutron stars 

CCOs Weakly-magnetized young neutron stars in centers of supernova remnants ~ 10° |Constant thermal-like soft X-ray emission. 

Pulsars Spinning, magnetized, young neutron stars 3320? |Pulsating radio emission, pulsating X-rays, gamma-ray emission. 

Magnetars | Very young neutron stars with extremely strong magnetic fields. Divided into soft gamma repeaters (SGRs)/|30 © SGRs: hard X-ray flares, giant bursts. AXPs: persistent and pul- 

and anomalous X-ray pulsars (AXPs). sating X-ray emission, transient radio emission. Additionally: fast 
radio bursts. 
Low-Mass X-ray Binary systems (M.S Mo) 
Bursters |Accreting neutron stars, which depict thermonuclear X-ray bursts 118¢  |Type-I thermonuclear X-ray bursts, burst oscillations. 
Recycled millisecond pulsars 
RMSPs _ |Rotation-powered millisecond pulsars. Like regular pulsars but spun up by angular momentum gain from the]472 ° —_|Millisecond-period radio pulsations, X-ray and gamma-ray 
accretion. emission. 
AMSP Accretion-powered millisecond pulsars. 18° X-ray pulsations. Sometimes also thermonuclear bursts. 
TrMSPs [Transitional millisecond pulsars. MSPs transitioning between accretion-powered and rotation-powered states|3 / X-ray pulsations, X-ray flaring, transient optical emission. 
Spider binaries 

Black RMSPs with irregular radio pulsations. Very low companion masses of Me < 0.06 Mo 38° X-ray emission, radio emission, optical measurements of the 

widows companion. 

Redbacks |RMSPs with irregular radio pulsations. Hydrogen-rich companion with Me 2 0.1 Mo 22€ X-ray emission, radio emission, optical measurements of the 
companion. 

Other Binary systems 

HMXBs |Accreting neutron stars with high-mass (Me Z 10 Mo) companions [100] [constant X-ray emission, X-ray pulsations 

Merging |Coalescing binary systems composed of a neutron star and another compact object (NS-WD, NS-NS, or|/4 £ precursor emission (?), gravitational waves, short gamma-ray 

compact |NS-BH). These systems can also be observed when they coalesce and merge. bursts, kilonova (seen in IR, optical, UV). 

binaries 

Environments 

PWNs Pulsar wind nebulae. The region inside a supernova remnant powered by a wind from a central pulsar. 33 + 28 "|Broadband emission from radio to X-rays and gamma-rays. Cos- 
mic rays. 

SNRs Supernova remnants. Large-scale structures that are the result of supernova explosions. The region is bounded] 294 ' |Multi-band electromagnetic emission, cosmic rays. 


by an expanding shock wave. 


Notes: 


“http://www.iasf-milano.inaf.it/~deluca/cco/main.htm (list of observed CCOs); 
’https://www.atn£.csiro.au/research/pulsar/psrcat/ (pulsars catalog; [91 ]); 
“http://www.physics.mcgill.ca/~pulsar/magnetar/main.html (list of observed magnetars); 
4https://personal.sron.nl/~jeanz/bursterlist.html (list of observed Type-I bursters); 
*“https://blacksidus.com/millisecond-pulsar-—catalogue/ (millisecond pulsar catalog); 


Í Patruno et al. 2022 [112]; 


€ GW170817, GW190425, GW200105, GW200115; 
h PWN + PWN w/o coinciding pulsar http: //www.physics.mcgill.ca/~pulsar/pwncat .html (PWN catalog; [77]); 
‘https://www.mrao.cam.ac.uk/surveys/snrs/snrs.data.html (list of SNRs; [54]); 
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Laboratories of gravitation 


The extreme density of the matter that is concentrated inside neutron stars makes them ideal laboratories of 
gravitation at the extreme limit because they can strongly deform the spacetime metric [98, 141, 114]: the 
high compactness of the star will curve the spacetime, whereas the rapid rotation can twist it. These effects 
create complex spacetime metric deformations, which modify the visual appearance of neutron stars. De- 
tailed modeling of these effects allows us to estimate the radius and mass of neutron stars from observations. 

Another frontier in the gravitational physics of neutron stars is the gravitational wave observations from 
merging neutron star binaries [14]. Modeling these systems relies on understanding and exploiting the dy- 
namical, time-dependent side of gravity: as the coalescing, space-time-deforming neutron stars merge, they 
will simultaneously vibrate spacetime itself. The resulting gravitational waves can then be detected and com- 
pared to general relativistic hydrodynamic simulations to constrain the properties of the neutron star matter. 
In practice, this allows measuring the tidal deformability of the ultra-dense matter of the coalescing neutron 
stars. 


Space-time deformations 


The high density of the neutron star leads to strong space-time distortions. The exterior space-time metric of 
a static, non-rotating spherical object is given by the Schwarzschild metric [98]’, 


dr? 


ds? = —(1—u)dt? + i +7°d0* +r’ sin? 0dọ?, (1) 


—u 
where ¢ is the coordinate time, r is the radial coordinate (defined so that the area at a time fixed time is 47r), 
O is the latitudinal angle, @ is the longitudinal angle, and u is the compactness parameter, 


Rs _ 2GMns 
Rys  Rysc? ’ 


(2) 


u= 


G is the gravitational constant, and Rs ~ 2.95(M/Mo) km is the Schwarzschild radius. A neutron star 
is so compact that the corresponding Schwarzschild radius, Rs ~ 4.1km, is slightly smaller than the star’s 
physical radius, Rys /Rs ~ 3. Hence, the star is on the brink of collapsing into a black hole. 

The neutron star compactness parameter, u, measures the severity of space-time distortions. Its value is 
somewhere between u ~ 0.3 — 0.5, leading to strong dilation effects of the coordinate time tf (« y 1 —u 0.8) 
and stretching of the radial coordinate (« 1/ vI =u 7 1.2) of the space-time close to the star’s surface. 


7 We follow the Misner et al 98 sign convention and use a metric signature of (—, +, +, +). 
8 Sometimes it can be more convenient to use the isotropic Schwarzschild metric of the form 


2 


1-4 iN’ 
ds? (=) dt? 4 (1 5) (dr? +F (d0? + sin? 0d¢?)), 6) 
2 


where F is the isotropic radial coordinate and # = M/F is the isotropic compactness. Surfaces of constant time in this metric 
are conformally flat; therefore, all the angles are represented without distortion. However, in this case, the measured distances 
depend on the location. See [105]. 
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Escape velocity from on top of such an object is very large, Vese = ,/2GMns/Rns ~ 1.93 x 10!°cms7!, and 
is comparable to the speed of light, ~ c/2. 


Rotating stars 


Rotation of the compact star, and hence the large angular momentum, will lead to additional shear-like 
distortions of the space-time metric [114]. In this case, we need to specify the exact latitude where the effect 
is studied; a typical choice is an equator with the corresponding equatorial radius Re, defined such that the 
circumference of a star (using the Schwarzschild metric radial coordinate) is 477R,, as measured in the local 
static frame. The radius at the pole is smaller, Rp < Re. 

For a star rotating with an angular velocity Q, the moment of inertia is Ins ~ 10 gcm*. Angular mo- 
mentum is then Z = IQ. The dimensionless angular velocity is 


= Q Q?2R3 
Q= = e (4) 
Qk GMys 


where the Newtonian mass shedding limit (Kepler limit) is Qk = y GM /R3. Observed neutron stars have 
Q ~ 0.01 — 0.1 resulting in strong rotational spacetime effects. 

Rotation can be incorporated into the space-time metric by expanding around the small parameter Q. The 
first-order expansion corresponds formally to the Kerr-like metric, where an extra angular velocity term of 
a local inertial frame appears. This is known as the frame-dragging effect, where the spacetime is rotating 
together with the star (i.e., in practice, velocity v + (v — vfa), where vfa is the frame-dragging velocity). 
Second-order expansion is formally related to Hartle-Thorne-like metric [60], which, in addition to the 
frame-dragging effects, has (two) quadrupole moments entering into the metric description.’ [25, 27, 17, 3, 
105] 

The rotation will also modify the actual shape of the star, making it oblate— the rotating star develops a 
bulge on the equator. The shape is numerically found to be close to an oblate spheroid, [99, 3, 105, 152], 


Q 2 
R(0) = Re b- (3) (0.788 — 2.06) cos” 6 | , (5) 


where R(@) is the radius at a colatitude 6 (angle from the pole towards the equator). 


Radiation from the star’s surface 


The strong stretching of the metric in the radial direction modifies the apparent visual look of neutron 
stars for distant observers. The most prominent effect is the gravitational light bending around the star that 
increases the apparent radius (that we would measure by looking at the star from a distance) to be 


° These are the multipole moments of the energy density and pressure. Their value depends on the coordinate system, whereas 
a single coordinate-invariant quadrupole moment can be obtained as a combination of the two [113]. 


Fundamental physics with neutron stars 11 


oo o1 02 03 04 05 06 07 00 02 04 06 08 10 12 14 
F/Fo F/Fo 


Fig. 1 Example images of rapidly rotating neutron stars. The colormap shows the radiation flux, F œ (1+ )~+, from rota- 
tionally distorted neutron stars as a function of the non-red-shifted, non-boosted flux, Fy (z = 0). Left panel shows an image 
of a typical massive neutron star (Mys = 2Mo, Rys = 13km) spinning with a frequency f = 600Hz being observed with an 
inclination angle (w.r.t. to the spin axis) of i = 60°. Right panel shows the visual appearance of the same neutron star with 
an extreme spin frequency of f = 1400Hz close to the mass-shedding limit. The ray-tracing is performed with the general 
relativistic ray-tracing code Arcmancer, [123]. 


Ro = i (6) 
Similarly, the dilation of the time component modifies the energy of the emanating photons. This effect, 
gravitational redshift, is also connected to u, and is given as 


1 
1 


l+z= (7) 


= 


It is a measure of how much, for example, the frequency changes, v¢/Vv; = (1 +z)! (where Vy is the 
observed and v; initial photon frequency), as a photon climbs out from the gravitational well created by the 
neutron star. The typical effect is a redshift, v; > vf, where the frequency is reduced during the climb. 

The photon energy is also modified by the Doppler-like boosting factor, 


_ vi-# 
o> 7 geal’ (8) 


where 


RQ 
0) = —— sin 0, 9 
BQ) = = ©) 
is the surface velocity of the rotating star, and ¢ is the angle between the photon momentum and the sur- 
face velocity vector, p = B(—sin@,cos(@),0) at the emission point. The Doppler boosting of the radiation 
can lead to a blue shifting of the radiation, vr > v;, causing an increase in the frequency, as photons are 
“slingshot” away from the surface like being thrown out from a rapidly rotating cartwheel. 
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The total frequency (or energy) change is the product of the gravitational and Doppler-boosting factors, 
vf/vi = 6/(1 +z). Observed flux scales as Fg œ (v;/v;)? and energy-integrated flux as F œ (v¢/v;)*. The 
emergent flux is visualized in Fig. | for two different example cases. Finally, we note that all the radiative 
processes should be defined in the co-rotating frame of the star (i.e., for an observer located at the surface). 

To summarize, the emission calculations typically take into account space-time deformation effects up to 
second-order in rotation, which are composed of: 


i gravitational potential (and hence gravitational redshift), 
ii Doppler-boosting, which originates from the rotating surface, 
iii change in the local emission angle (affecting the limb-darkening and brightening effects), 
iv time transformations (because the integration of finite-sized areas needs to take place simultaneously), 
and 
v oblate shape of the star. 


To a lesser extent, we also need to take into account the following: 


vi quadrupole-deformations of the spacetime (modifying the gravitational redshift term), 
vii frame-dragging effect (reducing the angular velocity close to the equator), and 
viii quadrupole deformations which modify the differential emission area element (dS = R*dcos @d@ for a 
flat spacetime metric). 


Typical neutron star images are constructed using the so-called S+D approximation (Schwarzschild met- 
ric and Doppler boosting) [115, 128, 127], where the observer’s polar coordinates are connected to the 
star’s co-rotating spherical coordinates (algebraically) and quantities computed at the surface of the star are 
Lorentz-boosted into the static non-rotating observer’s frame. The method can be refined by taking into 
account the oblate shape of the star [27, 99, 85, 105, 136, 24, 86]. 


Pulse profile modeling 


If the surface of the neutron star has a spot that is located off the rotational axis (in comparison to the 
whole surface just emitting uniformly), it produces periodic pulsations with the spin frequency of the star. 
Information in these pulse profiles can be used to constrain the geometry of the rotating system together with 
the mass and radius of the neutron star. 

Creation of such spots requires some physical effect that breaks the axisymmetry (along the azimuthal 
direction) on the surface and makes an azimuthal patch of the surface brighter (or, in theory, dimmer, cor- 
responding to cold spots that produce flux deficit as a function of the phase). Three mechanisms for spot 
creation are typically considered: heated magnetic polar caps in isolated pulsars, accretion-heated spots in 
binary systems, and dynamic thermonuclear burning regions. 

Heated magnetic polar caps of isolated pulsars. In isolated pulsars, the magnetospheric activity above 
the surface leads to acceleration of charged particles that are launched towards the surface [e.g., 135]. The 
in-falling “rain” of charged particles leads to heating the star’s surface- a hot spot. The opposite magnetic 
pole of the star can house a second antipodal spot. There are observational implications that some of such 
polar caps might not be circular but resemble more arc-like structures (banana-shaped spots) [131,96]. Such 
complex structures can be created, for example, if the magnetic field of the neutron star is not a dipole but 
has higher-order multipole components due to a more tangled internal field. 

Accretion-heated spots of neutron stars in (low-mass X-ray) binary systems. Accretion-heated spots 
are generally labeled as originating from accreting millisecond pulsars. In these systems, the surface is heated 
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by the material from the companion star, which falls into the neutron star. The spot is formed as the magnetic 
field of the neutron star channels the material into a stream that hits the surface and heats it. 

Thermonuclear burning fronts. The thermonuclear burning of accumulated material on top of accreting 
neutron stars (in low-mass X-ray binary systems) can sustain a coherent burning region, similar to hurricanes 
on Earth’s weather layer. The dynamics of these storms are still not very well understood: The burning fronts 
can lead to spots that drift (i.e., the spot is slowly moving on the surface) and expand (i.e., the burning front 
propagates into the fresh fuel), making their modeling very complicated. 


Irrespective of the origin of the spots, the information in the observed pulsating emission — pulse profile 
— can be used to measure the radius and mass of the neutron star housing the spot. The simplest assumption 
for an emitting spot is a circular emission region at a co-latitude 0; with a half-angular size p (measured 
in degrees using great-circle distances). During one rotation of the star, the spot first becomes visible to 
the observer, then crosses the visible surface, and finally disappears behind the star. This leads to a phase- 
resolved pulse profile where the detected emission increases, peaks, and decays. Depending on the observer’s 
inclination angle i and on the spot’s colatitude, the spot can disappear behind the star, experience a partial 
occultation, or always remain visible. 

As an example, NASA’s NICER (Neutron Star Interior Composition Explorer) mission has used the 
pulse profile analysis to constrain the radius and mass of the rotation-powered millisecond pulsar PSR 
JO030+0451 to be R ~ 12.7+1.2km and M ~ 1.3440.16Mo, as measured by [131]; or, R ~ 13.0t13km 
and M ~ 1.44+0.15Mo, as measured by [96]. These NICER observations also indicated that emission ar- 
eas were, in fact, more consistent with three elongated arcs instead of just two antipodal circular spots. This 
might imply that the magnetic field structure is more complex than the regularly assumed dipole field. Simi- 
larly, the NICER collaboration has constrained the radius of a massive pulsar PSR JO740+6620 with a known 
mass of M + 2.08 +0.07Mo to be R & 12.4413 km, as measured by [132]; or, R ~ 13.717$ km, as measured 
by [97]. For comparison, a similar pulse profile analysis was performed with RXTE (Rossi X-ray Timing 
Explorer) data to constrain the radius of the accretion-powered millisecond pulsar SAX J1808.4—3658 to be 
R ~ 11.9 +0.5km (assuming mass M = 1.7 Mọ) [136]. These measurements agree with neutron stars with 
a radius of R ~ 12km. 


Gravitational waves 


So far, we have discussed static gravitational effects like the gravitational bending of light rays or the redshift 
of photon energy. Recent pioneering observations of gravitational waves from neutron star mergers are a 
good example of dynamic gravitational phenomena where spacetime changes as a function of time [14]. 

A typical toy picture of a gravitational wave depicts the wave on a vibrating elastic membrane with 
a surface wave-like perturbation moving on it. This conforms to the geometric interpretation of gravity 
as a spacetime curvature, i.e., general relativity. Alternatively, we can imagine the gravitational waves as 
a consequence of special relativity in which a change in the gravitational field takes time to be felt at a 
distance; any perturbation (even of spacetime) will propagate with a maximum velocity of c. 

Imagine a general field with strength decreasing as an inverse square of the distance, « r~?. Such a field 
can be visualized by a familiar picture of radial field lines originating from the source of the field. If such 
a source is uniformly moving, the field lines will still have the same radial form. However, if the source is 
accelerating, there will be a propagating “dividing zone” between a moving radial field line region close to 
the source (expanding with a maximum velocity of c) and a further-away exterior region still depicting a field 
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of a non-moving source. The discontinuous zone dividing the interior and exterior regions has a width cAt, 
where At is the duration of the acceleration. The discontinuity can be understood as a wave traveling away 
from the accelerating source. This forms an alternative basis for understanding the nature of gravitational 
waves. 

Mathematically, the simplest object that can launch such a pulse is a time-varying monopole, X = 
dx /dt # 0, where X is, for example, the electric charge (X — q, where q is an electric charge; relevant 
for theories of electromagnetism) or mass (X — M, where M is mass; appropriate for theories of gravi- 
tation). However, the conservation of electric charge and conservation of mass prevents the existence of 
time-varying monopoles. So any perturbation traveling on a field generated by any kind of “change of a 
monopole” is forbidden. 

Let us next consider dipoles, the second-most simple mathematical object in line. In electromagnetism, 
an electric dipole (i.e., a pair of balanced positive and negative charges) is given as P} = } qisi, where 5; 
is the distance of the charge from the center of the charge cloud. Acceleration (for example, a rotating 
motion) corresponds to a non-zero second-time derivative: È, = OP, /dt? =Yqia; Z 0, where a; £ 0 is the 
acceleration of the ith charge. As is well known, this leads to electromagnetic dipole radiation. We return to 
this form of energy loss in the electromagnetism section. 

Gravitation turns out to be more complex since momentum is always conserved (in a closed system), 
and so By = 0? Py /dt* = $ Mia; = 0. To make a system emit gravitational radiation, time-varying higher- 
order moments are required. The subsequent mathematical entity in the multipole expansion after the dipole 
moment is the quadrupole moment (e.g., four balanced charges arranged in the corners of a square). The 
quadrupole moment of mass distribution is defined as Iy = } M;s;®s5;, where ® denotes the tensor product of 
two vectors, resulting in a quadrupole moment being a rank-2 tensor that is described with a3 x 3 matrix. The 
corresponding physical quantity is the gravitational tidal field, g’ = Ag/Ad, which represents an observable 
relative acceleration (i.e., a force) between two displaced particles within a small displacement of Ad (i.e., 
the gradient of gravity). 

The dominant gravitational radiation originates from the fourth time-derivative of the quadrupole mo- 
ment, 7, which produces a change in the gravitational (tidal) field g per displacement s of g’ = 0g/ds 
(~ Ag/Ad, where 0g ~ Ag and s ~ Ad). This change is proportional to g’ ~ `T /r, where r is the distance 
from the source. Other, weaker modes scale as g' ~ I/r?, [/r*, Ï/r’, and T /r?; all falling off more rapidly 
than the dominant 7 /r term. The dominating mode contains terms of the form Msä, Mva, and Ma?; for os- 
cillatory motion, relevant for rotating systems with equal masses M moving over distance s with frequency 
f, all of these scales as œ Mf*s*. The amplitude of the gravitational radiation then scales as 


GM f*s? s 1 
GEN ~ GM 1 
2 ó r G Atr’ (10) 


where A is the wavelength. 
In gravitational wave physics, the most commonly used parameter characterizing the amplitude of the 
wave is not g’ (i.e., change in gravity per displacement) but the dimensionless strain, 


n=2f fear ~2S al (11) 


(i.e., change in displacement per displacement). The two integrals in the strain formula represent an instan- 
taneous change in displacement as a function of time. Therefore, h is twice the fractional change in the 
displacement between two nearby masses due to the gravitational waves. This can be understood geomet- 
rically as follows. Similar to electromagnetic radiation, this displacement is on a plane transverse to the 
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direction of the motion of the waves; i.e., the wave is a shear wave. It stretches the spacetime along one axis 
and squeezes the other axis orthogonal to it. The net distortion is twice as large as what the stretching or 
squeezing would independently give. Note also that h itself is not observable; only the second derivatives 
and higher of h (that produce acceleration) are detectable. 

The final scaling of the strain is 


h~ == 


— y. (12) 


TNG 


In the most standard, already observed scenario of two in-spiraling neutron stars, the rotating binary system 
corresponds to a non-zero quadrupole moment that acts as a cyclic source of gravitational waves. A neutron 
star with a mass M ~ 2Mo, inspiral velocity v ~ c (just before the merger), at a distance r ~ 6 x 10% cm 
(~ 20 megaparsecs) would result in a peak dimensionless strain of 107%. This kind of fluctuation in the 
strain was indeed observed in the GW170817 event [1]. 

Other possible mechanisms for generating gravitational wave radiation from neutron stars are, e.g., break- 
ing of axisymmetry in rotating pulsars; here a small mountain on top of the neutron star (corresponding to 
some defect of the crust) would also lead to non-zero gravitational quadrupole moment, and hence emis- 
sion of gravitational waves. The perturbation amplitude from rotating mountains is much weaker, and these 
signals have not yet been detected. 


Interpretation of gravitational waveforms 


Latest observations of gravitational waves from merging double neutron star binary systems offer a powerful 
way for constraining the properties of the ultra-dense matter inside the stars [14]. The gravitational wave 
emission observed from neutron stars is often split into the chirp and ringdown phases. The chirp phase 
occurs at the beginning of the inspiral as the two stars begin to coalesce toward each other but are not yet 
touching. In this case, the gravitational wave frequency is just twice the orbital frequency of the binary, 
tow © 2fx, where fx is the Keplerian orbital frequency. The contact frequency, the frequency reached at the 
contact between the two objects, is well-approximated as 


3 33/2 
=a (13) 
2G M 
where č = GM /Rc? is the average compactness, assuming identical stars with Mı = M3 = M and R; = 
Rz = R. The chirp frequency is observable for tens of seconds and peaks at faw ~ 1 kHz. During the chirp 


phase, the equation of state can be probed by measuring the dimensionless quadrupole tidal deformability 


(polarizability coefficient), 
2 REN’ 
= 3h (a) pa 


where kz is the quadrupole Love number. Its measured posterior distribution is currently consistent with 
A = 0 (corresponding to two black holes merging). However, since the system is a known double neutron 
star binary (and, in addition, because a gamma-ray burst was detected), we are confident that the GW170817 
event measured the tidal deformability of the ultra-dense matter inside neutron stars. This measurement rules 
out the softest equations of states that require large A > 0. 

As the two stars merge, they will either collapse directly into a black hole or remain, for a while, as a 
“joint” neutron star. The resulting star can be either a stable neutron star supported by the standard pressure 
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of the matter or a quasi-stable supra massive/hyper massive neutron star, supported (at least partially) by the 
centrifugal acceleration from the rotation of the object. Supramassive neutron star remnants assume a uni- 
form rotation profile, Q (r) œ r (where r is the radius from the center of the remnant), whereas hypermassive 
remnants have differentially rotating interiors. 

The cores of the two merging stars continue to rotate inside a common envelope for a while, resulting in 
a ring-down gravitational wave emission. The ringdown phase is short (the system loses angular momentum 
rapidly due to efficient dissipation caused by the dense surrounding cloud) and emits gravitational waves at 
high frequencies. The ringdown phase has three characteristic frequencies, Q1 < Q2 < Q3, with a distinct 
physical origin. The physical origin of the frequencies can be understood by a toy model describing a rotation 
of two balls connected by a spring. The balls mimic the cores of the initial neutron stars, which rotate in a 
common envelope. Such a system can have three characteristic configurations if the system is rotating and 
the spring is vibrating. When the balls are at their maximum separation, the system has the largest moment 
of inertia and, conversely, the smallest rotation frequency, Q;. At their closest separation, the system has 
the smallest moment of inertia and, conversely, the highest rotation frequency, 423. If there is no dissipation, 
the angular frequency of the system will oscillate between these two states. Time spent at a given frequency 
is ~ Q/Q, and so most time is spent at Q ~ 0, i.e., at the two extreme states, resulting in power-spectrum 
peaks at frequencies corresponding to Q; and 3. If the spring is dissipative, the system asymptotes to an 
average frequency, Q2 ~ (Q1 + Q3)/2. This leads to a time-growing third peak at a frequency of Qo. 

Realistic merger simulations confirm that the initial few oscillations of the ringdown show a gravitational 
wave spectrum with two frequency peaks at f = Q,/m and Q3/7. Realistic merger remnant configurations 
are highly dissipative, so they quickly develop a strong peak at f ~ Q2/z, with weaker sidebands at fı and 
fs. Especially the fı frequency correlates with the compactness, u = 2GM /Rc?, enabling a more accurate 
equation of state constraints in the future when observations of the ringdown phase become possible. 

Finally, we note that the (short) gamma-ray bursts occurring shortly after the neutron star mergers are also 
an important new avenue of research [47]. They are sources of electromagnetic radiation (so-called kilonova) 
and possible neutrinos. They are also speculated to be responsible for the heavy element production in our 
Universe since the merger fuses extremely heavy elements, and the gamma-ray burst expels them to the 
surrounding interstellar medium. 
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Laboratories of nuclear physics 


The density inside neutron stars is so high, and the particles are packed together so tightly that quantum 
effects play an essential role in controlling the behavior of the matter [59]. This makes neutron stars exciting 
laboratories of nuclear physics. Furthermore, since the interiors of neutron stars are the densest form of 
observable matter known in the Universe, they provide a pathway into new studies exploring particle physics 
and quantum chromodynamics — something previously thought possible only with particle colliders. This 
line of research of the neutron star interiors is encapsulated into the unknown equation of state of the matter, 
P(p,T,...), describing a pressure P as a function of thermodynamic properties of the matter like density, p, 
temperature, T, and so on. 


Dense matter inside compact objects 


When the matter is compressed from ordinary matter towards the ultra-dense states inside neutron stars, 
different forces play a role hierarchically. First, the pressure is supported by the thermal motions of the 
gas. When compressed further, the charge repulsion between negatively charged electrons and positively 
charged protons becomes important. Both of these effects are explainable by classical physics and are used 
to describe the matter and equation of state inside regular stars. 

For compact objects, like neutron stars and white dwarfs, quantum mechanics is involved when explaining 
the behavior of matter. First, for dense enough matter (such as those found inside white dwarfs or outer 
layers of neutron stars), the electrons start to be packed so densely that quantum mechanics influences their 
motions via the so-called degeneracy pressure. When the compression is continued, the protons and neutrons 
will follow, providing a degeneracy pressure (such as in the core of neutron stars). The equation of state at 
these densities is still unknown. Beyond even the densities reachable in the interiors of neutron stars, we 
know that the final support of an ultra-dense matter is given by quarks via their degeneracy pressure (such 
conditions could perhaps exist inside the most massive neutron stars) [5, 6]. 


Degeneracy pressure 


Nuclear matter behaves in a somewhat analogous manner to liquids: the particles composing the liquid 
experience attraction between their neighbors but strongly resist compression when tightly packed. This 
resistance originates from zero-point motions of electrons (electron degeneracy pressure) and later on from 
zero-point motions and mutual interactions of neutrons (neutron degeneracy pressure). 

The degeneracy pressure is based on Pauli’s exclusion principle that forbids two identical fermions (spin- 
1 particles) from occupying the same quantum state. This introduces a quantum mechanical pressure as 
fermions repel their neighbors with similar quantum states. We can derive the resulting pressure from 
Heisenberg’s uncertainty principle, which asserts that position and momentum can not be simultaneously 
determined to be better than 


AxAyAzAp,ApyA p: = h°, (15) 


18 Joonas Nattila and Jari J.E. Kajava 


1089 
1038 
1037 4 
1036 


1035 


Mass M (Mo) 


1034 


Pressure P (dyne cm~?) 


1033 


1032 


1014 1018 106 9 10 11 12 13 14 15 16 
Density p (g cm~?) Radius R (km) 


Fig.2 Examples of the different equation of state models and the corresponding neutron star mass-radius relations. Each EOS 
has a one-to-one mapping between the density-pressure and mass-radius planes. 


where h is the Planck’s constant, and AX? Ap? is the 6D phase-space volume of the particles, where loca- 
tion vector x = (x,y,z), and momentum vector p = (px, Py, pz). This means that there exists a zero-point 
momentum, called Fermi-momentum, pr. 

The Fermi-momentum can be associated with a Fermi-energy, Ep. For non-relativistic case (p < mc) it 
is Ep = De /2m. Here we consider Fermi energy for a (cold) 3D space with N fermions filling the available 
energy states from the lowest up to the Fermi energy. It can be approximated similar to the energy in a 
square-well potential of size L, 


O nR 
~~ 2m 


2 
TE N E. (16) 
* » 2! Imi2 


Ep 


where 7, is the radius in n-space, and n; are the principle quantum numbers in 3D. The total energy is 
obtained by integrating all states inside the Fermi sphere, 


nhe Tn nH 


Eo = 276 a ò 4rr?dr, n= (17) 


Here the factor of 2 comes from the fact that we have two different spin states for fermions. 
Finally, the pressure is obtained from the standard thermodynamic relation (assuming a volume of V = 


L), m 
342 
Be OE tot _t he (3N (18) 
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assuming constant entropy (S) and particle number (N). We can generalize the expression into a form [59], 


pD P, x3 (xr <1) (19) 
Yap Pt (1), 


where P, = mec? [90° Àg ~ 1.5 x 10°dyncm™? is a typical pressure, Ac =ħ/mec = 3.86 x 107!! cm is the 
(reduced) Compton wavelength, x, = pp /mec is the dimensionless Fermi-momentum, and Yap = 3 or Yap = 


5 is the polytropic index for non-relativistic (x, « 1) or ultra-relativistic (x, >> 1) gas. Note that x, « ni! 9 x 
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Fig. 3 Equation of state, P(p), of the dense matter inside neutron stars. The upper horizontal axis gives the same relation 
as a function of radius, P(Rszy,1.4), for a M = 1.4Mo star with the SLy equation of state. Different regions inside the star 
are highlighted, including: atmosphere, for which the equation of state is given by the ideal gas law (various temperatures of 
T = 10°, green curve; 5 x 10°, blue curve; and 107 K, black curve are shown), ocean, where the plasma is Coulomb liquid, outer 


crust, where the equation of state is split between non-relativistic degenerate electron gas (P œ p5/3; magenta dashed line) and 


relativistic degenerate electron gas (P œ p“ 3; dashed blue line), and inner crust and the core where the P(p) is still uncertain 
(various models are shown with different colored curves). Location of the neutron drip, pyp (where neutrons start to leak out 
from the atomic cores) and saturation density, Pn (corresponding roughly to densities inside nucleons), are also shown with red 
dotted vertical lines. 


p 1/3. The transition between non-relativistic and ultra-relativistic regimes occurs at x; ~ 1 corresponding to 
p ~10°gcm-?. 

This theoretical treatment provides a pathway to understanding the nuclear physics of the cores of neutron 
stars. It is easy to see that degenerate neutron gas will behave similarly to degenerate electron gas since we 
can replace the mass me —> mp. Repulsive nuclear interaction, however, also contributes to the pressure at 
very high densities, rendering the equation of state in the neutron star core, p > 10'4gcm~?, unknown. 
Multiple competing nuclear equation of state models exist at these densities (see Fig. 2). Structure of a 
neutron star with one of these models, SLy (Skyrme Lyon model; [44]), is demonstrated in the previous 
Fig. 3. 

Most models assume fully nucleonic material (i.e., interactions between the neutrons). Some models also 
assume that the material is enriched with hyperons, nucleons containing strange quarks or that the stars are 
composed entirely of strange-quark matter (mixture of up and down and strange quarks), making them, in 
fact, strange stars, in comparison to neutron stars. The strange star models produce very soft equations of 
states which are generally disfavored by current observations of the heaviest neutron stars of M ~ 2Mo [e.g., 
43]. 

In summary, the state of ultra-dense matter is still open at the densities corresponding to the interiors 
of neutron stars. Eventually, for high enough densities, the material will decompose into quarks, with an 
equation of state given by the perturbative quantum chromodynamics. The equation of state in different 
density intervals is then composed of the following: 
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i ideal gas law, (P = NKT /V , where kT is the gas temperature) for p < 10° gcm~?, 


ii non-relativistic degenerate electron gas (P œ p>! 3) up top ~ 106 gem’, 


iii relativistic degenerate electron gas (P œ p*/ 3) up top ~ 10!! gem~?, 
iv degenerate neutron gas with additional, unknown, nuclear interactions, and, eventually, 
v quark matter with pressure given by quantum chromodynamics, at p = 40p, (where p, is the saturation 


density). 


The equation of state is also visualized in Fig. 3. Especially the unknown equation of the state of the matter 
inside the core is under active research. 

Observations of neutron stars allow an unrivaled opportunity to put definite constraints on the state of 
the ultra-dense matter. Various methods exist for constraining the mass, both mass and radius, and in some 
cases, the state of the matter indirectly. For recent reviews, see [109, 150, 95, 81] and especially [41] that 
focuses specifically on the various methods for mass and radius measurements. We have already discussed 
the possibility of using the pulse profiles (and their dependency on the gravitational potential) to provide 
constraints for the equation of state. Other promising methods include constraining the size of the emission 
directly from the X-ray emission emanating from the surface. Various methods used to constrain the equation 
of state of ultra-dense matter inside neutron stars are listed and described in Table 2. 


Thermal-like emission from isolated neutron stars 


The most important part of interpreting the electromagnetic emission from neutron stars is understanding the 
physics of the atmosphere that re-processes all the radiation going through it. The atmosphere of a neutron 
star is a thin, ~ 1 to 10cm height, a slab of hot plasma. To model it, equations of radiative transfer are 
solved, and the atmosphere is (typically) assumed to be in a radiative and hydrostatic equilibrium. This kind 
of modeling bears a great similarity to modeling regular stellar atmospheres. 

The atmosphere is thought to be composed of light elements. It is typically assumed to be in a local 
thermodynamical equilibrium (LTE) so that the equation of state follows an ideal gas law, P = nkT and Saha 
equation can be used to calculate the ionization levels. Heavier elements are expected to sink to the bottom 
due to the large gravitational potential (i.e., the atmosphere has a very short sedimentation timescale); hence 
the atmosphere is often taken to be composed of the lightest elements only, such as hydrogen or helium. For 
accreting systems, the composition is often taken to coincide with the elements in the companion star’s outer 
envelope. 

The most direct application of the atmosphere models is the interpretation of the thermal-like emission 
from central compact objects (CCOs), point-like X-ray sources in supernova remnants, using the atmo- 
sphere model spectra. These CCOs are young (t < 10*yr) neutron stars with modest X-ray luminosities of 
Ly ~ 10°? erg s~!. They have thermal spectra consistent with temperatures of kT ~ 0.1 to 0.5keV. Interest- 
ingly, recent observations indicate that some CCOs are better explained by the carbon composition of the 
atmosphere. Using carbon atmosphere models, Ho & Heineken [66] constrained the radius of a neutron star 
in Cas A to be between R ~ 10— 14km. Similar measurement was done for HESS J1731-347 by Klochkov 
et al. [149, 78] who constrained the mass and radius to be M = 1.55 + 0.3 Mọ and R = 12.4193 km, respec- 
tively. 


Method 


Radio pulsars in binary systems 


Table 2 Methods for constraining the neutron star equation of state 


Description 
Mass measurements 


Accurate timing of radio pulsars allow the determination of orbital parameters (orbital period Por, 
inclination angle sini, the mass of the companion M,, and total mass of binary Mr = Mys + Me) for 
neutron stars in binary systems. For nearly edge-on systems, the propagation of the signal through 
the gravitational potential of the companion allows breaking the degeneracy in the mass (Shapiro 
delay). 


Notable measurements 


PSR J161-2230: M = 1.908 + 0.016Mo [43, 48, 
11], PSR J0348+0432: M = 2.01+0.03M. [7], 
PSR J0740+6620: M = 2.08 +0.07Mo [49]. 


Optical measurements of binary 
systems 


Thermonuclear X-ray bursts 


Optical measurements of a companion star rotating around a neutron star allow the determination 
of orbital parameters (P5;p, sini, Me, and Mr = Mys + Me) of the system. Additional information 
via i) Shapiro delay, ii) observations of spectral lines of the companion star, or iii) modeling of 
eclipses (black widow systems) can be used to break the degeneracy between the mass and incli- 
nation. 

Simultaneous mass and radius measurements 


Time-resolved X-ray observations of cooling of neutron star atmosphere after thermonuclear X-ray 
bursts allow determination of R, M, distance D, hydrogen mass fraction X, and other system pa- 
rameters by comparing the cooling track to theoretical model calculations. Distance measurements 
help narrow the constraints. 


Black widow system PSR J0952-0607 M = 2.35 + 
0.17Mo [133]. 


Burster 4U 1702—429 R = 12.4+0.4km and M = 
1.9+0.3Mo [107]. 


X-ray pulse profile modeling 


Modeling of phase-folded X-ray pulse profiles allow constraining the parameters of the gravita- 
tional potential (compactness GM /Rc?). Higher-order effects like Doppler boosting and relativistic 
aberration of angles allow breaking of the degeneracy between R, M, and observer inclination an- 
gle sini and spot co-latitude cos 65. 


PSR J0030+0451 R © 12.7 £ 1.2km and M © 1.34 
0.16Mo [131] or R ~ 13.01} km and M ~ 1.44 + 
0.15Moọ [96]. PSR J0740+6620 with mass M ~ 
2.08 + 0.07Mo and R ~ 12.4+låkm [132] or R % 
13.7+7$ km [97]. 


Thermally-emitting qLMXBs 


Modeling of X-ray spectral observations of quiescent (non-accreting) X-ray binary systems al- 
low constraining the size of the emitting region (~ R? /D*) and relativistic effects carry additional 
information about compactness which allows weak constraints on the mass M too. Distance mea- 
surements help narrow the constraints. 


Statistical combination of five sources yielded R = 
9.1713 km for M = 14M, (90% confidence) [57]. 
Similar analysis of 8 sources yielded R = 10 — 14km 
for M = 1.4Mo [145]. 


Thermally-emitting isolated 


neutron stars 


Constraints from fast spin 


Modeling of X-ray spectra observations of isolated neutron stars allow constraining the size of the 
emitting region (« R? /D7), and relativistic effects carry additional information about the compact- 
ness, which allows weak constraints on the mass M too. Distance measurements help narrow the 
constraints. 

Other constraints 


Measurements of neutron stars with fast spin allow constraining the maximum possible values of 
M and R by requiring the star not to break up (mass shedding limit.) 


CCO HESS J1731—347 R = 12.4793km and M = 
1.55 +0.3Mo [78]. 


PSR J1748-2446a rotating f = 716.356 Hz [62] limits 
R < 16km (for M = 2Mo). 


Tidal deformability measure- 
ments 


Modeling of gravitational wave signals from mergers carry information on the tidal deformability 
of the material that can be translated to M and R constraints. 


GW170817 allowed constraining the tidal deforma- 
bility [1], which in turn translates to R < 14km for 
M= 1.4 Mo. 
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Thermonuclear X-ray bursts 


Thermonuclear (type-I) X-ray bursts are a phenomenon unique to neutron stars (for a recent review, see [50]). 
They offer a more complex alternative for measuring the neutron star radius by comparing the time-resolved 
X-ray spectra to theoretical atmosphere models. In addition, thermonuclear burning on top of neutron stars 
is a rich phenomenon with many unknown physical features. 

The bursts are triggered by unstable thermonuclear burning of matter accumulated onto the neutron star 
envelope from a binary companion. During the burning, hydrogen and helium fuse into heavier elements 
through the hot CNO cycle and triple alpha reactions, sometimes triggering the rp-process (rapid proton cap- 
tures) that can generate elements up to Tellurium [139]. Enormous amounts of energy are released quickly, 
particularly for helium-rich fuel: triple alpha reactions are not limited by beta decays that limit the CNO 
hydrogen burning. The most energetic bursts can have rise times that are less than an ms [70], suggesting 
that on rare occasions, the neutron star envelope can also detonate rather than deflagrate (i.e., the burning 
front can expand both super- and sub-sonically), and cause a relativistic outflow that drives the outer neutron 
star envelope into the interstellar space. 

Observations of X-ray bursts probe several key fundamental physics problems. In bursts where the rp- 
process is important, the (often poorly known) nuclear reaction rates of various proton-rich isotopes play a 
key role in the morphology of the X-ray light curves [39]. This has the potential to make the thermonuclear 
bursts also laboratories of nuclear physics, rivaling any Earthly nuclear burning physics experiment on nu- 
clear reaction rates. The extreme parameter regime also enables studies of more exotic proton-rich isotopes 
that can be present during the burning. 

The burning ashes generated during the X-ray burst — in addition to the possible, stable nuclear burning 
[138] — settle down to the neutron star ocean. At the bottom of the ocean, chemical separation occurs, 
wherein heavier elements crystallize into a solid crust, and lighter elements such as oxygen and carbon 
preferably remain in the liquid ocean [68]. This transitional layer of the star is where so-called super-bursts 
are triggered, which are powered by nuclear burning of carbon [37, 38, 146]. In the accreted crust [82], 
the chemical composition is further altered through pycnonuclear reactions that also heat the crust [58]. 
However, the neutron-rich nuclei in the outer crust undergo repeated cycles of electron captures and beta 
decay. The resulting neutrino losses effectively thermally decouple the ocean/atmosphere from the hotter 
inner crust [140]. Also, very recently, a new type of nuclear burning regime has been proposed, the so- 
called “hyper bursts”, that are triggered possibly by unstable burning of oxygen or neon at densities of p ~ 
10!! gcm~? near the neutron drip [1 11]. It is practically impossible to obtain direct observational constraints 
deeper down from these densities. Therefore, to probe the properties of the inner crust and the ultra-dense 
neutron star core, we must place constraints on the macroscopic properties of the star: their masses and the 
radii. 

We can place constraints on the neutron star masses and radii by analyzing quasi-thermal light emitted 
by their photospheres. During X-ray bursts, the entire atmosphere heats up, such that the photospheric tem- 
peratures can reach almost 3 keV, and then it rapidly cools down as the fusion reactions start waning down 
at the burning depths. The most energetic bursts reach the Eddington limit, which results in photospheric 
expansion, and they act as standard distance candles [79]. Moreover, these photospheric radius expansion 
bursts are expected to launch radiatively driven winds [1 10, 56], which may enrich the interstellar medium 
with the burning ashes [159, 69, 74]. 

Multiple techniques have been used to constrain neutron star masses and radii. The majority of them — 
such as the touchdown method [40, 108] or the cooling tail methods [148, 151, 106] — rely on the fact 
that their X-ray spectra can be adequately modeled with a simple black body model. Another approach is 
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modeling X-ray bursts using sophisticated atmosphere models, whose properties depend on the neutron star 
mass and radius, among many other variables [107]. 

All of these techniques, however, have relied upon the assumption that the environment where the X-ray 
bursts occur remains unaltered by the explosion of light and mechanical energy exerted by the burst-driven 
winds. Recent observational data indicate clearly that this is not the case the majority of times; the X-ray 
bursts alter the properties of the accretion discs surrounding the neutron stars [161, 162, 164, 75, 42]. The 
accretion flow in the disc conversely seem to influence the spectral properties of the bursts [72, 73]. The 
burst-disc interactions severely limit the use of X-ray bursts in making precise mass-radius measurements. 
Typically, their impact can only be inferred from studying large numbers of bursts rather than seeing clear 
imprints in individual bursts. In addition, the atmosphere models that have been used so far ignore many 
important factors, such as rapid stellar rotation that can cause the star to become oblate and modify the 
energy distribution of the radiation via strong Doppler boosting. 

But the potential of constraining dense matter equations of state models using X-ray bursts is enormous. 
For example, [107] constrained the neutron star’s radius in one X-ray burster to within 400 meters. This target 
was, however, ideal in many ways: the mass accretion rate was uncharacteristically low, and the source was 
in the hard spectral state, facilitating optimal measurement conditions for the radius. Unfortunately, most of 
the other observed bursts/bursters are less optimal in their behavior. 
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Laboratories of Electrodynamics 


Understanding the physics of neutron star magnetospheres requires the usage of both classical electromag- 
netism and quantum electrodynamics (QED) [9, 10, 29, 23]. From the perspective of electromagnetism, the 
neutron star is just a rotating magnetized ball in a vacuum. This kind of simplification, however, only gets 
us so far because the strong magnetic field also means that electromagnetic radiation (photons) and plasma 
particles start to a couple in the framework of QED. 


Spindown power of magnetized balls 


What can we expect from a dense ball of plasma (R ~ 10°cm) that is rotating rapidly (P ~ 1s) and spinning 
a huge magnetic field (B ~ 10!? G) around? The rotational energy of such a star is given by 


1 
E = IQP, (20) 
2 
where the moment of inertia J = 2MR? ~ 10* gcm?. Since isolated pulsars are observed to slow down, the 
corresponding spin-down power is 
E=12Q ~ 10” erg s™!. (21) 


for an angular velocity of Q =22/P ~ 6s7!, and spin-down rate of È ~ 107}. 

The simplest physical model for this slow-down is an electromagnetic torque acting on a rotating dipole. 
The angular momentum is carried away from the system by electromagnetic radiation, Poynting flux. The 
corresponding magnetic dipole energy loss rate is 

2(m,)? 2 B?Q*R® 


B=; a M sin? y ~ 10°°BY,Py terg s~}, (22) 


where m, = msiny, m = BR? is the magnetic moment, B the star’s polar magnetic field, and y is the 
magnetic dipole inclination angle with respect to the spin axis. The period and angular velocity are related 
simply as P = 27/ Q. In the last expression, the quantities are scaled as Q, = 10*Q. 

Equating the spin-down power and the magnetic dipole losses allows us to estimate the magnetic field. 


We obtain, 
[P | È 
~ 1012 


that agrees with our previous estimate of the strength of the magnetic field for characteristic pulsar parame- 
ters. 


Charges in the magnetosphere 


The above estimates demonstrate that the basic arguments based on electrodynamics lead to quantitative 
agreement. The standard pulsar model, therefore, consists of a rotating magnetized ball with a dipole field. 
Let us now refine this model to include charge carriers — plasma. 
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From the perspective of our “rotating ball model’, the magnetic field is spun by an excellent conductor. 
This rotation induces electric fields in the conductor. We can model the pulsar as a Faraday disk dynamo to 
estimate the order of magnitude of these electric fields. In this picture, we have a disk of radius R, rotating 
with an angular velocity Q. The disk is penetrated by magnetic fields of strength B ~ 10!? G. The rotating 
conductor experiences the moving magnetic field as an electric field, 


v 
E| =|- -xB| ~ —, (24) 
c c 


where v = Rf x Q is the tangential velocity of the outer disk rim. The corresponding voltage induced by 
the rotating disk (for a Faraday disk generator) is V ~ 10!5V. There is also a large current flowing in the 
“circuit”, I = P/V ~ 10'4 A (assuming spin-down power of P ~ 10°% erg s~!). 

The large electric field can lift charges from the star’s surface. Goldreich & Julian [53] showed that this 
induces a (minimum) charge density of 


V-E Q.B 
4n ~ 2mc` 


Pol = (25) 
These charges are pulled from the surface because the electromagnetic forces far exceed gravity. The cor- 
responding number density of plasma is ngy = Pcs /e ~ 10!2cm~?. An actual number density can be even 
higher, n = ngj, where .@ > 1 is the so-called multiplicity parameter. 

The charge density in the magnetosphere leads to a formation of a (minimum) conduction current of 


J =nev, (26) 


where ne = pp is the charge density required to screen the parallel electric 
E-B 
Ej = —— +0; (27) 
B 


in the opposite case, the electric field would lift the charges until it would become screened. 
The minimum Goldreich-Julian charge density [53] also has higher order relativistic corrections so that 
Poi —> Pas+ additional relativistic terms [102]. They modify the Goldreich-Julian charge density, 


(Q—@)B 


2 
2ac a 


Pal ~ — 
where @ ~ TAO is the angular velocity of the spacetime drag at a distance r from a rotating body. These 
corrections are, in reality, important to fully explain the plasma feeding into the magnetosphere. 


Force-free and magnetohydrodynamic solutions 


Much of the further progress in understanding the structure of the magnetospheres originates from numerical 
studies. This is understandable because the related equations are nonlinear and often require numerical tools 
to probe the solutions. The subsequent structure of the magnetosphere depends on how dense the plasma is, 
that is, if æ > 1. 
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Fig. 4 Schematic view of the pulsar magnetosphere. The neutron star in the center has an inclined dipole magnetic field w.r.t. 
to the spin axis of ¥ = 2/12. Open and closed magnetic field lines are shown with solid black curves; magnetic field lines 
forming the Y-point are depicted with dotted curves. The light cylinder is shown with blue dashed curves. 


Goldreich & Julian [53] introduced the idea of a force-free magnetosphere where the parallel electric 
fields are always screened. Here, the electromagnetic energy density is so large that all inertial, pressure, and 
dissipative forces can be neglected. This means that the pulsar can be considered immersed in a massless 
conducting fluid instead of a vacuum. Similarly, Michel [93] introduced a magnetohydrodynamical counter- 
part; this also assumes a large density, compared to ng). 

In both systems, the plasma starts to rotate rigidly together with the star. This rigid rotation becomes 
impossible at a sufficiently large radius. This is known as the light cylinder, where the magnetic field moves 
with the speed of light, and the plasma can no longer co-rotate. Therefore, the light cylinder radius is 

c 
Ric =J (29) 
It is Ric ~ 10° to 10!°cm (i.e. 10°R,) for P ~ 0.1 — 1s. Inside this region, r < Ryc, the magnetic field is 
primarily dipolar (higher order moments decrease faster with r than dipole moment). Outside of it, r > Ric, 
all the field lines are open, and the flux is conserved so that Br? = const.; then, B œ r7!. 

The problem of magnetospheric structure can be formulated in a type of Grad-Shafranov equation that en- 
capsulates the aligned rotator; the equation is familiar from magnetic confinement [16] with Grad-Shafranov 
equation being its relativistic counterpart [137, 94]. Contopolous et al. [35] solved this equation by iterative 
technique assuming E -B = 0 and E? — B? < 0 (implying subluminal E x B drift velocity). The important out- 
come of this exercise was that the last closed field at the equator has a Y-type neutral point. This means that 
the return current flows in a current sheet at the boundary of the closed zone where the magnetic field flips 
sign. The result has since then been produced by others confirming its universal nature [55, 92, 156, 144] 
The corresponding spindown energy losses of the aligned rotator were numerically confirmed to be 


Qu? 
3.0 


E=k (30) 


C 


with k = 1 £0.1. 
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Evolving magnetic topology 


Important progress came when Spitkovsky [144] numerically solved the magnetospheric structure with time- 
dependent force-free simulations. This confirmed the presence of the current sheet at the equator and allowed 
the solution with an arbitrary (magnetic) inclination angle 7 > 0 (angle between the magnetic moment and 
rotation axis). Then, the generalized spindown is 

204 


Q 
z (1+sin? x) G1) 


É =k 


with k = 1 +40.1. This expression introduces the E œ (1 + sin? X) dependency on the inclination angle, %. 
The magnetic topology of an oblique rotator, % > 0, is just a rotationally distorted version of the simpler 
aligned rotator geometry, 7 = 0. Note that the solution still assumes a magnetosphere filled with the plasma 
of density sufficient to short out the parallel electric fields. 

This generalized spindown power formula suggests an increase in the losses with increasing 7. Therefore, 
the inclination angle should decrease with time for old pulsars. Since Ric increases as the pulsar spin downs, 
there must be a net conversion of open field lines to closed field lines on the spin-down timescale. This 
topological change requires the violation of the ideal MHD. This will then necessarily lead to reconnection 
[89]. The most likely location is at the Y-point and in the current sheet. The importance of reconnection 
was later confirmed with fully kinetic simulations that do not require us to impose the assumption of dense 
plasma background [119, 120, 34, 122, 121, 30, 31]. 


Mysterious Pulsar Radio Emission 


The spin-down luminosity, ~ 10°% erg s~!, leaks out from the pulsar magnetosphere and can act as free 
energy to produce observable emission. However, only a tiny fraction, 10~° to 1074 of that energy is known 
to end up in the radio band to generate the actual pulsar emission [see | 18, for a recent review]. Quoting 
J. Arons, this makes pulsars “dogs that don’t bark” [9]. The origin of this emission is still unknown, but 
it is thought to be connected to the pair avalanches (see next section) regulating the electric circuit of the 
magnetosphere. 

The radio emission is broadband, extending to as low as tens of MHz and up to 150 GHz. The energy 
output peaks at a frequency of a few hundred MHz with a specific energy flux of a few Jansky [118]. !° The 
emission is elliptically polarized with a high degree of polarization, reaching ~ 100%. The linearly polarized 
component dominates the emission. Circular polarization degree can reach up to ~ 10%. The polarization 
angle of many pulsars has a canonical swing, where the angle evolves as a function of the pulsar phase, 
tracing a (tilted) S-shaped trajectory. This angle sweep takes place independently of frequency, pointing to 
geometric effects that are well explained with a rotating vector model [130, 126]. 

Observations averaging over multiple pulses of radio emission have given the pulsars their status as the 
most accurate clocks in the Universe. Individual pulses, however, vary significantly from pulse to pulse, 
and many show high sub-pulse variability. Some pulsars have individual pulses that drift in phase, others 
sometimes switch off (nulling), and others abruptly flip between various states (moding). The origin and 
nature of these sub-pulse variabilities are still unknown. 


10 Jansky is a unit of energy flux used primarily in radio astronomy. It is 10~73 erg s~! cm~? Hz7!. 
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Pulsar wind nebulae 


So where does the rest of the large spin-down energy go? Most of that energy is carried away by electromag- 
netic Poynting flux and by the charged-particle outflow streaming out from the pulsar. But all that energy 
and charged particles are not just freely flowing into a vacuum; the neutron stars are surrounded by gas from 
their old stellar remnants. The interaction of that electromagnetic energy and charged particles generates a 
hazy, glowing turbulent blob of gas around the neutron star known as a pulsar wind nebula. This makes the 
nebulae huge calorimeters in the sky that trap, reprocess, and weigh the energy output of pulsars. 

As an example, the Crab nebula is probably the most well-studied nebula (and neutron star) of all time. 
Crab has a spin period of P ~ 33ms and a period derivative of P ~ 420 x 10~!>. Its age can be estimated 
from P/2P ~ 1000yr (dynamical age), matching well with the historical supernova in 1054AD. The nebula 
emits in a broad range of spectrum from radio to optical to X-rays to gamma-rays; injected powered of E = 
5 x 108 erg s~! is required to be channeled into the surrounding plasma to explain the observed luminosity of 
the nebula. These estimates are well in-line with the spin-down power of a magnetized ball with B ~ 10!*G, 
corresponding to a power of ~ 10°% erg s~!. 

To better understand the interaction of energy and the charged particles, the radiation and its generation, 
we need an understanding also of plasma physics. This is the focus of our next section. 
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Laboratories of plasma physics 


The magnetospheres of neutron stars are composed of charged particles — plasma. As the magnetospheric 
plasma cools down, it radiates its energy away and populates the medium with photons. The nonlinear 
dynamics of this plasma and radiation are governed by radiative plasma physics. In this section, we introduce 
the basic concepts of such systems. 


Standard quantum electrodynamic interactions 


The magnetospheres of regular radio pulsars are composed of a large magnetic field of B, ~ 10!*G. Rapid 
rotation of this field leads to a generation of large electric fields, as well; these electric fields can open 
large potential gaps in the magnetosphere. The resulting voltage accelerates particles to relativistic energies, 
and these high-energy particles produce luminous radiation. The interaction of the charged particles and 
the radiation field requires the coupling of radiative physics and electrodynamics. Therefore, the pulsar 
magnetospheres are prime laboratories of quantum electrodynamics [7 1, 22]. 

In practice, quantum electrodynamics describes the interaction of photons (electromagnetic radiation) 
and charged particles. The simplest two-body interactions between photons and electrons are:!! 


i pair creation (photons interacting with photons to create pairs), y+ Y — e7 +e”; 


ii pair annihilation (electrons interacting with positrons to create photons), e~ +e* > Y+Y; 

iii Compton scattering (photons interacting with pairs to exchange energy), e~ + y > et + y; and 

iv Coulomb scattering (pairs interacting with pairs to exchange energy), e* +e* — e™ +e*™ and e* + 
et et + eF", 


Magnetic field adds an extra flavor to these (lowest-order) interactions because it can be thought of as con- 
sisting of virtual photons. This enables single-body interactions (or, alternatively, two-body interactions 
composed of “real” against “virtual” species) such as 


v synchrotron radiation (pairs upscattering virtual photons), e+ + [B] + e* + y+ [B]; and 
vi single-photon (magnetic) pair creation (photons pair-creating with virtual photons), y+ [B] + e* + [B]; 


All of these processes are important in the neutron star magnetospheres. 


Pair cascades 


The magnetosphere of pulsars is prone to avalanches of pair creation [147]. They evolve as follows [23]. The 
twisting electric field (from the rotation of the star) induces a longitudinal electric field to the magnetosphere, 


E~ (QR/c)B. Any (primary) particle in this region is accelerated to ultra-relativistic energies, € >> 1, 


where £; = E,/ mec? is the particle’s energy, Ee, in the units of electron-rest mass meC2. 
The strong magnetic field of the environment, however, renders the synchrotron cooling time to be very 


short, 


11 We give a schematic description of these interactions by denoting electrons with e~, positrons with e+, and photons with 
y. Pairs are succinctly given as e*. Altered (or often energized) end-products are denoted with a prime, x — x’. The passive 
magnetic field (or virtual photons) is expressed with [B]. 
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~ 107s, (32) 


where @g = eB/mec is the particle’s gyro-frequency and re = e? /mec? classical electron radius. Since the 
synchrotron emission power depends on the perpendicular component of the particle’s velocity, Psynch % 
€2 sina, where sing is the angle between the particle’s momenta and the magnetic field. This forces the 
particles to move along the magnetic field lines; pairs in the magnetosphere can be imagined as beads sliding 
on a wire. 

The synchrotron cooling and the curvature of the magnetic field causes these electrons (or positrons) to 


emit gamma-ray photons, £y ~ 1, as they slide along the B field lines. Here €y = Ey/ mec? is the energy of 


. ‘ be cat 2 
the photon in the units of electron rest mass. The energy losses from the curvature radiation, P, = $ Se, 
C 


depend on the curvature radius, Re ~ Rys, and the energy of the electron, £e, and impose an upper bound on 
the particle’s energy. Balancing acceleration, eE I> and the curvature losses, P}, results in Emax ~ (R2E | /e) 1/4 


(i.€., €maxMmec? ~ 107MeV). Therefore, the magnetosphere is expected to be full of high-energy gamma rays. 

The gamma-rays, corresponding to energies €y ~ 1, single-photon pair-produce with the magnetic field, 
resulting in low-energy secondary pairs, €& ~ 1. The required magnetic field strength can be estimated by 
considering the energy gap between the adjacent Landau levels, A@g, and comparing it to the electron rest 


23 
mass energy, i@g = mec?. This results in Borp = mee ~ 4.4 x 10!3G. Before pair producing, the mean free 


12 


path of photons in the magnetosphere is 


1 ~ Re Ban 1 
Y 10 B & 


(33) 


This is ~ Re ~ R, for MeV photons. 

As a result, the neutron star magnetosphere is quickly filled with e™ that screen the initial voltage. The 
creation of secondary pairs is halted when the longitudinal electric field, Ej, is screened. Therefore, consid- 
ering a rotating magnetized ball in a vacuum is not enough since the magnetosphere is dynamic and capable 
of generating its own plasma. 

This nonlinear feedback mechanism causes cyclic bursts of pair creation where a single seed electron can 
trigger a generation of a huge number of pairs. The resulting charges are blown out relativistically into the 
surrounding space. After the magnetosphere has cleared from the pairs, the voltage starts to grow again. This 
leads to cyclic generation of bursts. 


Vacuum birefringence 


For extremely large magnetic fields, especially those of magnetars with B ~ 10!°G, even the properties of 
the surrounding vacuum itself become modified. This happens because of another quantum electrodynamic 
phenomenon called vacuum birefringence. 

Birefringence, in general, is a property of, for example, crystals in which the speed of light inside the 
crystal is dependent on the direction of the light with respect to the axis of the crystal. In crystals, this 
breaking of isotropicity can happen even if the crystal lattice is, for example, more tightly packed in one 
direction. This causes the crystal to have different refractive indices of light depending on the direction; 


12 Omitting a logarithmic function that weakly depends on the neutron star parameters. 
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birefringence means that instead of one regular, isotropic refractive index, nı, the medium has two indices, 
fi and n2. 

The vacuum around strongly magnetized neutron stars can become birefringent and develop a differing 
refractive index along and perpendicular to the local magnetic field direction, ny, and n | . The effect is caused 
by virtual electron-positron pairs of the vacuum background spontaneously appearing and annihilating. Like 
any charged particle, these virtual pairs also gyrate in the strong magnetic field. This creates an imaginary 
crystal lattice in the vacuum that modifies the propagation of electromagnetic waves through it. Most notably, 
the polarization of the radiation is forced to align with the ambient magnetic field. 

The strength of the effect can be quantified as [64], 


2 
OF B ae) 
nn = 307 G=) sin’ @, (34) 


where ær ~ 1/137 is the fine structure constants, and & is the angle between the direction of propagation 
and the external field. The magnetic field remains strong enough to a distance known as the polarization 


limiting radius, 
2/5 1/5 
B E 
7 Y E 
ra~ 10 (wz) (i) (sin œ)? cm. (35) 


This renders the radiation in the inner regions, r < rp, strongly birefringent and, therefore, also modifies 
the polarization of the radiation. If the inner field is coherent, like a dipole field, the ordered magnetic field 
means that the polarization direction of the radiation is also coherent. An observer integrating the polarized 
radiation over a large solid angle containing the face of the star will hence detect a large net polarization. 
Therefore, a high degree of polarization is expected for radiation from magnetars. 


Superfluid and superconducting interiors 


The material at the core of the neutron star is thought to be in a superfluid state [18]. The strong interaction 
force between the nucleons renders the charge-neutral neutrons into a superfluid state and the electrically 
charged protons into a superconducting state. The electrons are expected to remain in their regular state be- 
cause the corresponding superconducting transition temperature is very small compared to the temperatures 
relevant for neutron stars. 

The superconductivity can lead to drastic changes in the magnetic properties of the object. Most crucially, 
superconductors exhibit the Meissner effect, where magnetic flux is expelled from the superconducting re- 
gions. However, the diffusion timescale of flux repulsion from macroscopic regions for the core material is 
very long (tp ~ 10”s) owing to the enormous electrical conductivity (large density enables easy transport 
of charges). 

The protons in the core are thought to form a type-II superconductor. This state is characterized by a 
periodic array of quantized vortices of super-current, which is aligned parallel to the local magnetic field. !* 
The magnetic flux associated with an individual proton vortex is 


13 On the other hand, Type-I superconductors have fine-scale, alternative regions of normal material encapsulating magnetic 
flux and superconducting regions exhibiting the Meissner effect. The state of the superconductor depends on the ratio of proton 
coherence length to the penetration depth. 
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and the number of vortices per area is, therefore, B/@o, meaning that billions of super-current vortices can 
thread the core. For B ~ 10!?G, the vortex lattice has a lattice separation of roughly 5 x 10~!°cm. This is 
large compared to the inter-particle distance, meaning that the vortices are mesoscale in size. 

The neutrons in the core are also in a superfluid state. A salient feature of the superfluid state is that 
the angular momentum of the superfluid becomes quantized into vortices that carry the angular momentum. 
Therefore, the core is expected to be composed of multiple vortices carrying the angular momentum of the 
rotating star. These superfluid vortex lines are parallel to the rotation axis. Each neutron pair vortex supports 
h of angular momentum; multiple vortices are, therefore, required to support the total angular momentum. 

The existence of a superfluid component and the possibility of a superconducting core makes neutron 
stars the largest blobs of superfluid material in the Universe. They allow us to study superfluidity and super- 
conductivity at the extreme limit since the core can remain superfluid even up to temperatures of millions of 
Kelvin due to the immense pressure. 


Gliches and quakes 


The high accuracy of pulsar timing measurements makes it possible to measure even the tiniest change in 
the star’s rotation rate. These high-precision measurements have revealed a phenomenon called glitches. 
Glitches are a rare, short-duration pulsar timing phenomenon. They are seen as sudden positive jumps of 
rotational frequency, Av > 0 (determined as the difference between the frequency after and before the glitch), 
in the pulsar timing data, followed by a negative change of the slope, v < 0. The glitch magnitudes vary from 
Av/v~10~' to 1075 and they have been detected in over 100 rotation-powered pulsars and magnetars [46]. 

Most glitches are followed by a recovery phase, where the spin-down rate increases and tends toward the 
pre-glitch value. The recovery phase can be modeled with an exponential function, œ e™* /t with a typical 
relaxation time of ~ 100 days together with a longer time-scale component with a relaxation time-scale of 
~ 1000 days. 

Physically, the glitches are thought to occur because of a rapid transfer of angular momentum between the 
superfluid interior and the outer crust of the neutron star, driven by catastrophic unpinning of the superfluid 
vortices. [4, 61] The angular momentum mismatch is created by slowing down the star’s crust with electro- 
magnetic torques. The rotation of the superfluid interior of the star, on the other hand, remains unchanged 
because the quantized vortices remain pinned. After the rotation mismatch exceeds an unknown threshold 
value, the vortices are unpinned and begin to slide along the lattice plane. Migration of the vortices away 
from the axis of rotation results in a spin-down of the superfluid component, bringing the two components 
closer to co-rotation. 

Magnetars are observed to also anti-glitch, where sudden negative, Av < 0, rotation rate changes, i.e., 
spin-downs (rather than conventional spin-ups) of the star are observed [8]. These events are speculated to 
occur due to magnetospheric reconfiguration. 

Something needs to trigger the catastrophic unpinning of the superfluid vortices. The most plausible 
candidates are starquakes and catastrophic readjustments of the crust [134]. The quakes are expected to 
occur when the crust accumulates stress beyond a critical strain that it can sustain, leading to a mechanical 
failure of the crystal lattice [67]. Why or how such large stress can build up in the crust lattice is still 
unknown. 
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Seismic motions of the crust will also couple the pulsar interior and the magnetosphere. The quake is 
expected to launch an elastic shear wave into the crust lattice, and since the magnetic field is frozen in the 
crust and liquid core, it also leads to perturbation of the magnetosphere. The resulting crustal oscillations 
deform the magnetic field lines and excite Alfvén waves on the magnetosphere [e.g., 26]. These waves can 
also interact with each other and heat the magnetospheric plasma [e.g., 84, 104]. 


Giant bursts and fast radio bursts from magnetars 


A similar kind of star quake event, only more extreme, is thought to also power magnetar flaring phenomena. 
In this case, super-strong crustal quakes are proposed as triggers of giant bursts [45]. Similarly, smaller 
quakes could produce fast radio bursts as a side product of magnetospheric activity. 

Giant bursts. Giant bursts are massively energetic blasts of radiation. A typical giant burst consists of a 
short, ~ 0.2 sharp pulse of hard X-rays and soft gamma-rays of energies E ~ 100keV to MeV. The initial 
pulse is followed by a fainter, slowly-fading tail lasting minutes. In the case of the famous March 5th, 1979 
event, the tail also varied in intensity with a sinusoidal perturbation (with two differing peaks per cycle) and 
a period of 8s. These oscillations were seen for more than 20 cycles. 

The leading model that could power these flares is a large-scale magnetospheric reconfiguration — these 
models rely on the large magnetic energy reservoir of the magnetar and on the fact that energy stored in the 
magnetic field can be quickly released to explain the extremely short duration of the pulses. In this picture, a 
strong crustal quake launches Alfvén waves into the magnetosphere. The waves twist the field lines, causing 
them to inflate and expand. The inflated bundle over-twists and undergoes rapid magnetic reconnection. 
The magnetic reconnection leads to the creation of strong non-ideal electric fields, and turbulent plasma 
motions that can energize and heat the local plasma [e.g., 103]. The mechanism resembles the energization 
mechanism observed behind solar flares. 

The twisted, evolving bundle carries a strong electric current and a relativistic stream of particles along 
the arched magnetic line loop. The streaming charged particles quickly radiate their energy as synchrotron 
radiation and can furthermore Compton up-scatter the photons up to relativistic energies. These high-energy 
photons will pair-create, forming a trapped pair-plasma cloud. Cooling of this pair plasma in the strong 
magnetic field provides a mechanism for generating the ultra-luminous soft gamma-ray burst seen from 
the soft gamma-ray repeaters. Furthermore, the streaming plasma will also fall back to the star, heating the 
magnetic footprints. Self-consistent modeling of the giant bursts requires detailed plasma physics models in 
strong magnetic fields, and large radiation densities [2 1]. 

Fast radio bursts. Similar plasma-physics-related phenomena to giant bursts are the recently found fast 
radio bursts (FRBs) [36, 116, 117]. FRBs are short ~ Ims duration radio pulses with energies ~ 10° to 
10* erg s~! (translating to 50mJy to 100Jy) observed between radio frequencies from 400 Mhz to 8 GHz. 
The emission is strongly polarized, exhibiting linear polarization degrees up to 100%. FRBs are the biggest 
radio-emission-based enigma since the discovery of pulsars. The latest coincident observation of FRB, to- 
gether with a flaring magnetar, undeniably established a connection between the two (at least for some 
sub-class of FRBs). 

The extremely high brightness temperature required to explain the FRB emission requires a coherent 
electromagnetic emission mechanism. The high coherence means that many particles are required to emit in 
phase rendering the problem automatically into a plasma physics problem. Current leading coherent emission 
mechanisms are the synchrotron maser instability from magnetized collisionless shocks [143] and plasma 
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waves (magnetosonic fast modes) from collapsing, reconnecting current sheets beyond the light cylinder of 
the magnetars [87, 90]. These and other emission mechanisms are reviewed in [88]. 


Extreme particles: cosmic rays, neutrinos, and more 


Finally, it should be mentioned that neutron stars are thought to be the origin of many high-energy particle 
events like cosmic rays, neutrinos; and, more speculatively, even sources of dark-matter candidates like 
axions. Observations of all of these can be used to probe fundamental symmetries in physical laws because 
their energies are unreachable by Earthly laboratories. In addition, they help us better understand the nature 
of particle acceleration in our Universe. 

Cosmic rays. Cosmic rays are relativistic protons or heavier atomic nuclei moving with velocities near 
the speed of light. The flux of cosmic rays is constantly raining down on us on Earth. They are typically 
detected when they slam into the upper layers of Earth’s atmosphere and dissipate their energy in showers 
of secondary high-energy particles (like X-rays, protons, electrons, alpha particles, pions, muons, neutrinos, 
and neutrons). The secondary particles have been first observed decades ago and are now routinely mon- 
itored. In the most extreme case, a dim fluorescent light from a cosmic ray exciting atmospheric nitrogen 
can be detected if the night is dark and moonless. On a more regular basis, they are detected by space- 
borne observatories such Fermi Gamma-Ray Space Telescope (detecting the secondary gamma-rays from 
the showers) or ground-based observatories such as HESS or HAWCK (detecting Cerenkov radiation from 
the cosmic rays entering the Earth’s atmosphere). 

The cosmic ray energy spectrum peaks at ~ 10? eV and extends into a broad power-law tail. The power- 
law tail has two breaks: first one, (named “knee”) at ~ 10! eV and another (named “ankle”) at ~ 10!8 eV. 
The lowest energy particles, < 10! eV, are thought to originate mostly from our own Sun (due to solar flares 
accelerating particles), while the intermediate energies, 10!°eV < Ecr < 10!5 eV, are galactic cosmic rays, 
accelerated by shocks in supernova remnants [5 1, 2]—neutron stars in action! Other more distant sources of 
high-energy cosmic rays could include, for example, active galactic nuclei (supermassive black holes at the 
centers of galaxies). 

The most energetic cosmic rays, named ultra-high-energy cosmic rays (UHECRs), can have kinetic en- 
ergies up to 10% eV. In comparison, the LHC (Large Hadron Collider) on Earth can only probe energies up 
to ~ 10!7eV. These extreme particles are composed of the lightest elements, most likely pure protons or 
light nuclei like helium cores. However, each carries a kinetic energy equivalent to a baseball (m, ~ 150g) 
traveling with a velocity of vp ~ 100kmh7! ~ 2800cms™!, 


Ekin = Tmp? ~ 11508 x (2800ems~!)? ~ 6 x 108 erg ~ 4 x 10%eV. (37) 
The UHECR particle is extremely relativistic since its energy is Eucr = YUCRMpC’, so that Yucr ~ 10!° for 
Eucr ~ 10° erg. Their energy is so high that they are theoretically expected to quickly lose it by scattering 
against the low-energy photons of the cosmic microwave background (Greisen-Zatsepin-Kuzmin limit; GZK 
limit). This means that they cannot be produced by, for example, the most distant quasars or active galactic 
nuclei but must instead originate from somewhere nearby. The big question is, what mechanism can accel- 
erate protons to such extreme Lorentz factors in our neighborhood? This mystery is complemented by the 
latest observation of PeV photons from the Crab [83] that also strongly challenges our picture of the particle 
acceleration and its limits. Neutron stars are plausible engines for powering the particle acceleration in many 
cosmic ray generation models. 
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Neutrinos. Another class of extreme particle is the neutrinos — light-weight (my < 0.120eV) subatomic 
particles that interact only via weak interaction and gravity. Neutrinos come in three leptonic flavors: electron 
neutrinos (Ve), muon neutrinos (Vu), and tau neutrinos (v+). They oscillate between the different flavors; a 
neutrino with a specific flavor is composed of a quantum superposition of all three flavors and can randomly 
change its appearance via flavor oscillations. They are produced in radioactive decay, nuclear reactions, or 
as secondary particles by cosmic rays interacting with Earth’s atmosphere. 

Most of the neutrinos detected on Earth come from the Sun. However, some neutrinos are confirmed to 
have a cosmic origin. These “cosmic neutrinos” have been directly observed from a core-collapse supernova 
explosion SN 1987A [65], making supernovae and neutron stars prime candidates for neutrino sources. The 
primary mechanism for neutrino emission in the supernova explosions (i.e., non-degenerate matter with 
density p < 10’ gcm~*) is the annihilation of electron-positron pairs into neutrino-antineutrino pairs, e~ + 
e*t — v+¥ (where v is the neutrino and V is anti-neutrino) [19]. At high temperatures, T > 10°K, the 
volumetric energy loss rate from the pair annihilation is (see Ref. [101] for discussion), 
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Similar mechanisms could be at play in ultraluminous X-ray sources if they are interpreted as systems with 
super-Eddington accretion channeled onto magnetized neutron stars [15, 13]; the high mass flux, production 
of an accretion shock, and the resulting shock-heated high temperatures would render the accretion columns 
on top of these neutron stars into efficient sources of neutrinos — neutrino pulsars [101]! 

Detection of astrophysical neutrinos holds great promise for probing fundamental physics. In the event 
of a close-by galactic supernova, there is hope that we could accurately measure the neutrino energies and 
flavor spectrum. This would help us better constrain the core-collapse supernova physics and any high-energy 
non-standard interactions beyond the Standard Model of particle physics. 

Axions. Neutron stars are also speculated to work as probes of various dark matter candidate particles 
[157]. The most intriguing of these is the axion, an extremely light (ma < 1eV) elementary particle, pos- 
tulated to resolve the strong CP symmetry problem in quantum chromodynamics. Axions couple to the 
electromagnetic fields and can change the behavior of the electrodynamic processes in the magnetosphere. 
Most notable of these couplings is the alteration of the physics at the magnetospheric gap region, E - B Æ 0, 
which is associated with a strong voltage, efficient production of gamma-rays, and therefore also the possible 
origin of axions [129]. 


Summary 


Research of neutron stars encompasses many branches of modern physics, from gravitational physics to 
nuclear and particle physics, from electrodynamics to plasma physics. What makes neutron stars unique is 
that they exhibit extreme behavior in each one of these fields, often offering novel insight into the limits of 
our physical theories. 

We started by presenting a (short) history of neutron stars and their discovery. For a more thorough histor- 
ical discussion on neutron stars, we refer the reader to [163]. We then reviewed the basic arguments behind 
the studies of gravitation, nuclear & particle physics, electromagnetism and & quantum electrodynamics, and 
plasma physics in the context of neutron stars. Each of these fields is still full of mysteries and unexplained 
phenomena that need to be solved to understand the physics and astrophysics of neutron stars better. 
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In this chapter we have provided a quick precursor to the field of fundamental physics of neutron stars, 
and for the interested reader more thorough reviews of these topics can be found from the following studies. 
Excellent sources of more information for studies of gravity include [98, 114]. The physics of emergent 
radiation is discussed in [124] and the effects of light-bending in [105, 152]. Notably, a simple model of 
light bending in curved spacetime is developed in [20, 125]. Studies of the neutron star equation of states 
are reviewed in [109, 81]. Measurements of neutron star radii are discussed in [95, 150, 41]. The recent 
gravitational wave observations are presented, for example, in [14]. Studies of neutron stars and pulsar 
magnetospheres are reviewed in [29, 23]. Pair cascades and quantum electrodynamic effects are discussed 
in [154, 153]. Magnetars and their observations are reviewed in [76], and plasma physics of fast radio bursts 
are well-discussed in [88]. 


Cross-References 


e Low mass X-ray binaries (Bahramian & Degenaar) 
e Fundamental physics with black holes (Carcia) 

e Isolated neutron stars (Borghese) 

e Low magnetic field (Di Salvo, Burderi, Iaria) 

e High magnetic field (Mushtukov) 
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